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The dynamics of the Luttinger model after a quantum quench is studied. We compute in detail 
one and two-point correlation functions for two types of quenches: from a non-interacting to an 
interacting Luttinger model and vice-versa. In the former case, the non-interacting Fermi gas 
features in the momentum distribution and other correlation functions are destroyed as time evolves. 
In the infinite-time limit, equal-time correlations are power-laws but the critical exponents are found 
to differ from their equilibrium values. In all cases, we find that these correlations are well described 
by a generalized Gibbs ensemble [M. Rigol et al. Phys. Rev. Lett. 98, 050405 (2007)], which assigns 
a momentum dependent temperature to each eigenmode. 



I. INTRODUCTION 

Most of the theoretical effort in the field of strongly 
correlated quantum systems over the past few decades 
has focused on understanding the equilibrium properties 
of these fascinating systems. For instance, achieving a 
complete understanding of the phase diagram of rather 
"simple" models like the two-dimensional fermionic Hub- 
bard model still remains a huge challenge. Nevertheless, 
however important these endeavors are, understanding 
the phase diagram and the equilibrium properties of the 
phases of strongly correlated systems will not certainly 
exhaust the possibilities for finding new and surprising 
phenomena in these complex systems, especially out of 
equilibrium. 

In classical systems, the existence of steady states out 
of equilibrium is well known. Very often, however, the 
properties of such states have very little to do with the 
equilibrium properties of the systems where they occur. 
Moreover, also very often their existence cannot be in- 
ferred from any previous knowledge about the equilib- 
rium phase diagram: They are emergent phenomena. 

One good example of a classical non-equilibrium steady 
state is provided by the appearance of Rayleigh-Benard 
convection cells when a fluid layer is driven out of equi- 
librium by a temperature gradient. Indeed, it is known 
that dissipation plays an important role in the forma- 
tion of these classical non-equilibrium states. However, 
different from classical systems, dissipation in quantum 
systems causes decoherence, which usually destroys any 
interesting quantum interference effects. Thus, although 
one may also wonder if non-equilibrium steady states 
can also appear when quantum systems are driven out 
of equilibrium, the study of non-equilibrium phenomena 
in quantum many-body systems has been regarded, until 
very recently, as a subject of mostly academic interest. 
The reason for this may be decoherence due to coupling 
with the environment, which is always present in most 
experimental realizations of large quantum many-body 



systems, and which prevents the observation of coherent 
quantum evolution for long times. 

However, the recent availability of highly controlable 
systems of ultracold atomic gases has finally provided the 
largely lacking experimental motivation for the study of 
non-equilibrium phenomena, leading to an explosion of 
theoretical activitv. Ultracold atomic gases are es- 

pecially interesting because they are very weakly coupled 
to the environment, thus remaining fully quantum coher- 
ent for fairly long times (compared to the typical duration 
of an experiment). At the same time, it is relatively easy 
to measure the coherent evolution in time of observables 
such as the density or the momentum distribution. Thus, 
theorists can now begin to pose questions such as: As- 
suming that a many-body system is prepared in a given 
initial state that is not an eigenstate of the Hamiltonian, 
how will it evolve in time? And, more specifically, will it 
reach a stationary or quasi-stationary state? If so, what 
will be the properties of such a state? How much memory 
will the system retain of its initial conditions? 

From another point of view, the problem described in 
the previous paragraph can be formulated as the study 
of the response of a system to a sudden perturbation in 
which the Hamiltonian is changed over a time scale much 
shorter than any other characteristic time scales of the 
system. In what follows, we shall refer to this type of 
experiment as a quantum quench. Quantum quenches 
are also of particular interest to the 'quantum engineer- 
ing' program for ultracold atomic gases. @D] The reason is 
that, if we intend to use these highly tunable and control- 
lable systems as quantum simulators of models of many- 
body physics [such as the two-dimensional (2D) fermionic 
Hubbard model mentioned above] , it is utterly important 
to understand to what extent the final state of the quan- 
tum simulator depends on the state in which it was ini- 
tially prepared. In particular, one is interested in finding 
out whether the observables in the final state state can be 
obtained from a standard statistical ensemble (say, the 
microcanonical, or the canonical ensemble at an effective 
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temperature). If this is so, one would speak of thermal- 
ization. If this does not happen, then how much memory 
does the system retain about its initial state beyond the 
average energy E — {H)l 

We would like to emphasize that the above questions 
are not a merely academic. Indeed, ultracold atomic sys- 
tems allow for the study of non adiabatic dynamics when 
the system is driven between two quantum phases such 
as a superfluid and a Mott insulator. [HI Also, in a 
recent experiment, J43] it has been shown that a faithful 
realization of the Lieb-Liniger model l44j exhibits absence 
of thermalization. In other words, when prepared in an 
non-equilibrium state, the experimental system reached 
a steady state that cannot be described by any of the 
'standard' ensembles of Statistical Mechanics. This ab- 
sence of thermalization seems to be a consequence of the 
integrability of the Lieb-Liniger model, that is, the exis- 
tence of an infinite number of independently conserved 
quantities. This conclusion was backed by the theoreti- 
cal analysis of Rigol and coworkers, [TD] who noticed that 
the non-equilibrium dynamics of an integrable system is 
highly constrained. Thus, based on numerical simula- 
tions for the Tonks-Girardeau limit of the Lieb-Liniger 
model, these authors conjectured that the long-time val- 
ues of some observables should converge to those ob- 
tained from a generalized Gibbs ensemble, which can be 
constructed using the maximum (von Neuman-) entropy 
principle. [TUJ USJ |3S] The conjecture was first analyt- 
ically confirmed by analyzing an interaction quench in 
the Luttinger model by one of us. [5] Later, it has been 
also found true in other integrable models: Cardy and 
Calabresse studied a quench in a Harmonic chain, [15] 
Eckstein and KoUar analyzed the Falikov-Kimball model 
in infinite dimensions. |18| and the 1/r Hubbard model in 
one dimension. |T9j Moeckel and Kehrein ,20] studied an 
interaction quench in the Hubbard model in infinite di- 
mensions by a flow equation method, and found that the 
system reaches an intermediate non-thermal state. Fi- 
nally, recent numerical studies also have suggested that 
lack of thermalization may even persist in the absence 
of integrability in one-dimensional systems, [13] or that 
it may occur only certain parameter regimes of non- 
integrable models. [Hj 

However, it can be expected that, [26] for a rather 
general choice of the initial state, along with a situation 
where there are few conserved quantities, the system will 
lose memory of most of the details of the initial state and, 
after it reaches a steady state, the expectation of many 
experimentally accessible observables such as the particle 
density or the momentum distribution, will look essen- 
tially identical to those obtained from a standard thermal 
ensemble. [65] Indeed, this is what seems to be observed 
in the vast majority of the experiments with ultracold 
atomic gases. However, for experimental many-body sys- 
tems in general, it is hard to quantify whether this will 
be always the case. We should take into account that 
(except perhaps in the case of ultracold atomic gases) 
the exact form of the quantum Hamiltonian is frequently 



not known with accuracy. And even when it is known, it 
is not always possible to tell a priori whether the system 
is integrable or even if it has other conserved quanti- 
ties besides the ones assumed by the standard thermo- 
dynamic ensembles. As a possible experimental check, 
we can say that, provided the final result is largely inde- 
pendent of the particular details of the preparation of 
the initial state and that the observations agree with 
those obtained from a thermal ensemble, we can say that 
thermalization has occurred. Indeed, some recent nu- 
merical evidence. [23] supplemented by the extension to 
many-particle systems of a conjecture known as 'eigen- 
state thermalization hypothesis' (first introduced in the 
context of quantum chaos[47]), seems to indicate that 
lack of integrability will in general lead to thermalization 
(in the sense defined above). Indeed, Reimann (26J has 
recently analytically demonstrated that, under realistic 
experimental conditions, equilibration will be observed 
in an isolated system that has been initially prepared 
in an non-equilibrium mixed state. Nevertheless, even 
if the issue of thermalization for non-integrable systems 
may have been settled, other questions such as the de- 
tails of the transition from the integrable case (which 
thermalizes to a generalized Gibbs ensemble) to the non 
integrable case, which thermalizes to the standard micro- 
canonical or Gibbs ensemble for large enough systems, 
are questions that are still far from being completely un- 
derstood. [5S] 

In this article, we will not try to answer the diffi- 
cult questions posed in the previous paragraph. Instead, 
we focus on analyzing the quench dynamics of a rela- 
tively well-known one dimensional model, namely, the 
Luttinger model (LM). A brief account of the results in 
this article has been already published elsewhere. [S] In a 
future publication, we shall also deal with another closely 
related model, the sine-Gordon model. [48] 

The LM Hamiltonian can be represented as a quadratic 
form of creation and destruction operators. The long 
time behavior following a quench for Hamiltonians of 
this form has been recently considered by Barthel and 
Scholwock. These authors provided some general con- 
ditions for the appearance of dephasing and steady non- 
thermal states. This question has been also taken up 
recently by KoUar and Eckstein. [19] However, since the 
Luttinger model may be relevant to experiments using 
ultracold atomic gases (see Sect. IV I or numerical simu- 



lations, it is important to obtain analytical results. The 
simplicity of this model also allows us to test in detail a 
number of general results. [T^ H5] 

The rest of this article is organized as follows: In 
Sect. [ll| we discuss the evolution of some simple cor- 
relation functions of the Luttinger model. We consider 
the case where the interactions between the fermions is 
suddenly switched on, and the reverse situation, that 
is, when the interaction in suddenly switched off. In 
Sect. |III[ wc discuss how the infinite-time behavior of 
some of the correlation functions following a quantum 
quench can be obtained from a generalized Gibbs ensem- 
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ble. We also discuss some observables for which this con- 
jecture fails. The experimental relevance of our results is 
briefly discussed in Sect.|IV[ along with other conclusions 
of this work. Finally, the details of some of the lengthier 
calculations are provided in Appendixes A to C. 



II. THE LUTTINGER MODEL 

The Luttinger model (LM) describes a one-dimensional 
(ID) system of interacting fermions with linear disper- 
sion. It was first introduced by Luttinger [50] but its 
complete solution was only later obtained by Mattis and 
Lieb,[5T] who showed that the elementary excitations of 
the system are not fermionic quasi-particles. Instead, 
Mattis and Lieb introduced a set of bosonic fields de- 
scribing collective density modes (phonons) of the sys- 
tem, which are the true elementary low-energy excita- 
tions of the LM. The methods of Mattis and Lieb bear 
strong resemblance to the early work of Tomonaga [55] 
on the one-dimensional electron gas. Extending the work 
of Tomonag weel as that of Mattis and Lieb, Luther 
and Peschel [53] computed the one and two-particle cor- 
relation functions in equilibrium, thus showing that cor- 
relations exhibit (at zero temperature and long distances) 
a non-universal power-law behavior signaling the absence 
of long-range order. Later, Haldane [54H56] conjectured 
that these properties (i.e. collective elementary excita- 
tions exhausting the low-energy part of the spectrum as 
well as power-law correlations) are distinctive features of 
a large class of gapless interacting one-dimensional sys- 
tems that he termed (Tomonaga-) 'Luttinger liquids'. Us- 
ing the modern language of critical phenomena, the LM 
can be understood as a fixed point of the renormalization- 
group for a large class of gapless many-body systems in 
one dimension: the equilibrium properties at low temper- 
atures of many ID systems are universal in the sense that 
they can be accurately described by the LM. However, in 
this work we shall be concerned with non-equilibrium 
properties of the LM, and because the latter can in- 
volve highly excited states, we shall make no claim for 
universality. The precise conditions under which the 
results obtained here apply to real systems that are 
in the Tomonaga-Luttinger class should be investigated 
carefully in each particular instance (see discussion in 
Sect.[!v|. 

The Hamiltonian of the LM can be written as follows: 



p,a—r.l 



(1) 
(2) 

(3) 



^4 - ^ E 54(9) ■ Mq)Jo.{~q) : . (4) 

q.a—r,l 

Here (p) and ipl^ [p) are fermion creation and annihila- 



tion operators at some momentum p and J a (q) = X]p • 
V'q(p + q)4'a{p) '■■ The index a — r,l refers to the chi- 
rality of the fermion species, which can be either right 
(r) or left (/) moving; the symbol : . . . : stands for nor- 
mal ordering prescription for fermionic operators. This 
is needed to remove from the expectation values the infi- 
nite contributions arising from the fact that the ground 
state is a Dirac sea,i55] namely a state where all single- 
particle fermion levels with p < are occupied for both 
chiralities. This defines a stable ground state (at the non- 
interacting level), which in what follows will be denoted 
by |0). 

A. Bosonization solution of the LM 

In this section we briefiy review the solution of the 
LM. The Hamiltonian in Eqs. ([l]) to Q can be written 
as a quadratic Hamiltonian in terms of a set of bosonic 
operators. [51] First we note that the density operators 
Jaiq) = Yip ■ '^iiP + qHa{p) ■ obey the following com- 
mutation rules: 

[Ja{q), Jfsiq')] = (^1^^ Sq+q,fi5al3, (5) 

which can be transformed into the Heisenberg algebra of 
the familiar bosonic operators by introducing: 

b(.q) = imM-q) - ^i-q)Mq)] ■ (6) 

where ^^(q) is the step function. Note that the q = 
components (known as 'zero modes') require a separate 
treatment since Jq,(0) = is the deviation, relative to 
the ground state, in the number of fermions of chirality 
a = r,l. However, rather than working with Nr and Ni, 
it is convenient to introduce: 



N = Nr+Ni J:^Nr-Ni, 



(7) 



which, since N,. and TV; are integers, must obey the fol- 
lowing selection rule (—1)^ = (—1)'^ when the Fermi 
fields obey anti-periodic boundary conditions: ipa{x + 
L) = —ipaix) (L is the length of the system). Therefore 



■aQ\p\ isapx 



Tpaip), being 
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+1 and p = 2{n — ^)n/L, where n is an integer, and 
flo ^ 0+.[55] 

The Hamiltonian Hlm can be expressed in terms of 
the bosonic operators introduced in Eq. ([6]): 



Ho = J2hvF\q\bHq)b{q) + ^{N' 



J' 



(8) 



H2^1Y1 92iq)\q\ [b{q)bi-q) + bHq)bH-q)] 



q^O 



2L 



(9) 



H, = Y: %4(9)kl b\q)biq) + ^^(AT^ + J2y (10) 



q^O 
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Ignoring the zero mode part, the above Hamiltonian 
has the form of Eq. (Al), with the following identifica- 
tions: uo{q) = vfIqI, m{q,t) = 54(g)kl, and g{q,t) = 
.92 (9)!'? I, and it can be therefore be brought into diag- 
onal form by means of the canonical transformation of 
Eq. (A3). Hence, the Hamiltonian takes the form of 
Eq. 1X51 with uj{q) = v{q)\q\, being v{q) — {[vp + 
g4{q)r ~ [52(9)]^}^^^, and 9 7^ 0. As to the zero mode 
contribution: 



2L 



2L ' 



(11) 



where vn = WF+54(0)+g2(0) and vj = WF+54(0)-g2(0). 
This defines the equilibrium solution of the LM. In the 
following sections we shall be concerned with the quench 
dynamics of this model. 



B. Suddenly turning-on the interactions 

Although it is possible to solve the general quench 
problem between two interacting versions of the Lut- 
tinger model, we shall focus here on the cases where the 
interactions described by H2 and _ff4 (cf. Eqs. [S] and 
[4]) are suddenly switched on (this section), and switched 
off (next section). Thus, in this section, we shall as- 
sume that we have made the replacements 32,4(9) — ^ 
52,4(91 i) — .92,4(9)^(0 Eqs.(3[4|. The Hamiltonian 
at times i > is therefore the interacting LM. In 
other words, using the notation introduced in Sect. [Aj 
Hf = Hq + H2 + H4 = i?LM, whereas the initial Hamilto- 
nian (for t < 0) is Hi — Hq. However, we note that, since 
both zero modes, J and N, are conserved by Hi = Hq 
and Hf = TJlm, their dynamics factors out, and we shall 
assume henceforth that we work within the sector of the 
Hilbert space where J — N — (this sector contains the 
non-interacting ground state, |0)). Therefore, from now 
on, we shall omit Hzm in all discussions. 

As to the initial state, we shall consider that, within 
the spirit of the sudden approximation, at i = the 
system is prepared in a Boltzmann ensemble at a tem- 
perature T described by 



Po = Pit - 0) 



(12) 



where Zq = Tr e~^'^^ . We shall further assume that the 
contact with the reservoir is removed at t = 0, and that, 
after the quench, the system undergoes a unitary as it is 
in isolation from the rest of the universe. 



Eq. ( A7 ) in the Appendix describes the solution to the 



interaction quench in terms of the modes that annihi- 
late the initial ground state |0), the solution itself is not 
particularly illuminating. To gain some insight into the 
properties of the system following the quench, let us com- 
pute a few observables. Amongst them, we first turn our 
attention to the instantaneous momentum distribution, 
which is the the Fourier transform of the one-particle 
density matrix: 

C^Ax,t) = {e^"''^^i^lix)MOK'"''^%. (13) 



where (• • • )o means that the expectation value is taken 
over the ensemble described by po (cf. Eq. (12)). The 



time dependence of the operators is dictated by H{, 
as described in Section |^ Notice that, since in gen- 
eral [H{, Po] 7^ 0, time translation invariance is broken, 
and the above correlation function is explicitly time- 
dependent. 

The time evolution of ^a{x) can be obtained using the 
bosonization formula for the field operator: [SSI [S3 EZ] 



(14) 



being r]r, rji two Majorana operators (also known as Klein 
factors) obeying {ria,rii3} = 2Sap, which ensures the an- 
ticommutation of the left- and right-moving Fermi fields 
(recall that Sa = +1 for a = r and Sq = —1 for a = I). 
The bosonic fields: 

(j)a{x) = Saipoa + -J^^o, + *a(a;) + ^ai^), (15) 

where [Na, ^op] — i5a,p, and, in terms of Fourier modes, 
/2_-\ 1/2 

^Jx)^ lim V — e-''""/2e"°?^6(s„g). (16) 

The details of the calculations of C^^{x,t) have been 
relegated to the Appendix [B] In this section we will 
mainly describe the results. However, a number of re- 
marks about how the calculations were performed are 
in order before proceeding any further. We first note 
that interactions in the LM are assumed to be long 
ranged. [50j EU [55] This can be made explicit in the in- 
teraction couplings by writing 32,4(9) = 52,4(9-^0): where 
the length scale i?o ^ i is the interaction range. Thus, 
just like system size L plays the role of a cut-off for 'in- 
frared' (that is, long wave-length) divergences, the in- 
teraction range, i?o plays the role of an 'ultra-violet' 
cut-off that regulates the short-distance divergences of 
the model. The results given below were derived assum- 
ing a particular form of the interaction (or regularization 
scheme) where the Bogoliubov parameter (cf. Eq. A4 ) is 
chosen such that sinh2/3(q) = ^ e"'"'^"'/^. Furthermore, 
we replaced z;(g) by w = w(0). Indeed, these approxima- 
tions are fairly similar to the ones used to compute the 
time-dependent correlation functions in equilibrium, [53j 
given that the expressions that wc obtain for the out- 
equilibrium correlators are fairly similar to those of the 
equilibrium correlations. [53] This regularization scheme 
greatly simplifies the calculations while not altering in a 
significant way the asymptotic behavior of the correlators 
for distances much larger than i?n.[67] 

Returning to the one-particle density matrix (cf. 
Eq. 13), we note that it can be written as the product of 



two factors: 



C^^{x,t)^d°}{x)K{x,t), (17) 

where C'^\x) is the noninteracting one-particle density 
matrix, and thus hr{x,t) accounts for deviations due to 
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the interactions. Hence, this factorization allows us to 
obtain the instantaneous momentum distribution func- 
tion as the convolution: 



fiP,t) 



dk 
2^ 



f^'>Hp-k)hrik,t), (18) 



where = {e^'^'P/'^ + is the Fermi-Dirac dis- 

tribution, and 



hr{k, t) 



-f-oo 



(19) 



Before presenting the results for the expression of the 
one-particle density matrix as well as the momentum dis- 
tribution at finite temperatures, it is worth considering 
the much simpler looking zero-temperature expression. 
We first discuss finite-size effects. For a system of size L 
we obtain: [S] 



Rn 



d{x\L) 

d{x - 2vt\L)d{x + 2vt\L) r'/^ 



d{2vt\L)d{-2vt\L) 



(20) 



where d{z\L) = L| sin(7r2:/L)| is the cord function and 

clp^ixlL) = i{2Lsm[TTix + iaa)/L]y^ (aq ^ 0+) the 
noninteracting one-particle density matrix. Notice that 
this result is valid only asymptotically, that is, for 
d{x\L),d{x ± 2vt) ^ Rq. Thus, we see that Ctp^{x,t\L) 
is a periodic function of time with period equal to tq = 
L/2v. This is in agreement with the general expectation 
that correlations in finite-size systems exhibit time recur- 
rences because the energy spectrum is discrete. Although 
the recurrence time generally depends on the details of 
the energy spectrum, in the LM the spectrum is linear 
Lu{q) ~ v\q\ and therefore the energy spacing between 
(non-degenerate) many-body states is Aq — 2'Khv/L. 
Hence, the recurrence time tq ^ 27r/Ao follows (the ex- 
tra factor of i is explained by the so-called light-cone 



effect, see further below). 



The recurrent behavior ex- 
implies 



hibited by the one-particle density matrix ( 13 
that, after the quench, the system does not reach a time- 
independent stationary state as time grows. A similar 
conclusion is reached by analyzing other correlations, 
such as e.g. the finite-size version of the density cor- 
relation function, [5] 



Cj^{x,t\L) = {e'"^"^Jr{x)Jr{Q)e-'"^"'^)^ 
(l+72)/47r2 



[d{x\L)Y 



[d{x~2vt\L)f [d{x + 2vt\L)f' 

(21) 



where 



Jr(x) 



\ ^ e 



^ ^ Jr{q) 1pl{x)ijr{x) : = ^l9^0r(x), 
9 



(22) 



is the density (also referred to as 'current') operator in 
real space. 

In the thermodynamic limit, L — > oo, the recurrence 
time Tq — L/2v — >■ +oo, and the system does indeed 
reach a time-independent steady state. In this limit, 
d(x\L) — >■ and the single particle density matrix 
becomes: [9] 



C^Ax,t> 0) = Gi"Hx) 



Rq 




x2 - {2vtf 


7^2 


X 




{2vtY 





Hence, 



hr{x, t) 



Rq 
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x^ - {2vtf 


X 




{2vtY 



7V2 



(23) 



(24) 



In order to understand the evolution of the momentum 
distribution, without actually having to compute it, it 
useful to consider the various limits of the above expres- 
sion, Eq. (23). First of all, for short times such that 
2vt <^ \x\, the function hr{x,t) is asymptotically just a 



time-dependent factor, [3] 



hr{x,t) ~ Z{t) = 



2vt J 



(25) 



which can be interpreted as a time-dependent 'Landau 
quasi-particle' weight in an effective (time-dependent) 
Fermi-liquid description of the system. In other words, as 
time evolves after the quench, we could imagine that the 
quasi-particle weight at the Fermi level is reduced from 
its initial value, Z{t = 0) = 1, to Z{t > 0) < 1. At zero 
temperature, this time-dependent renormalization of the 
quasi-particle weight reflects itself in a reduction of the 
discontinuity of the momentum distribution /(p, t) at the 
Fermi level (which is located at p = in our convention) . 
Therefore, at any finite time, the system behaves as if it 
was a Fermi liquid and therefore it keeps memory of the 
initial state, that is, a non-interacting Fermi gas. 
Yet, for t +00, hr{x,t) becomes a power-law: |5] 



lim hj.{x, t) 

t^oo 



Rq 



(26) 



and also does the momentum distribution. The behavior 
of the momentum distribution at different times is de- 
picted in Fig.[T] Interestingly, this time-dependent reduc- 
tion of the quasi-particle weight after being quenched into 
the interacting state has been also found in Ref . [20l when 
studying an interaction quench in the Hubbard model in 
the limit of infinite dimensions. In this case, however, 
the discontinuity remains finite even for t — >■ -|-c», which 
is different from the behavior of the LM, which is known 
to be a non-Fermi liquid at equilibrium. Some of these 
non-fermi-liquid features also persist in the quench dy- 
namics. 
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Z=l 



(a) 



Z(t) 



(b) 




FIG. 1; Schematic time evolution of the momentum distribu- 
tion f{p, t) at zero temperature, (a) At t = 0, the momentum 
distribution is that of non-interacting fermions, with a dis- 
continuity at the Fermi level (p = 0) Z = I. (b) At f > the 
discontinuity is reduced in a power-law fashion Z{t) ~ t~'' . 
(c) For t — > -l-oo the discontinuity disappears and the mo- 
mentum distribution exhibits a power-law singularity close to 
the Fermi level p = 0, f{p,t — >■ +oo) — ^ — const, x \p\'^ . 
However, the exponent characterizing the singularity is not 
the equilibrium exponent. 



facts: i) Because of the variational theorem, the ini- 
tial state {i.e. the ground state of the non-interacting 
Hamiltonian = Hq), is a complicated excited state 
of the final Hamiltonian H[ — Hi^m- H) Both and 
Hi are critical {i.e. scale free, apart from the cut-off 
Rq), thus, the system is likely to remain critical. To 
end the present discussion, we note that, in the litera- 
ture on Tomonaga-Luttinger liquids, it is customary to 



introduce the dimensionless parameter K 



-2,3(0) 



terms of which 7^ = {K — K~^)'^/4, which needs to be 
compared with the equilibrium exponent resulting from 
7cq = {K + — 2)/2. We shall use the parameter K 
in other expressions below. 

It is worth emphasizing that the particular evolution 
of the asymptotic correlations from Fermi liquid-like at 
short times to non-Fermi liquid-like at infinite time ex- 
hibited by the one-particle density matrix, is also found 
in other correlation functions. However, the idea that the 
system 'looks like' an interacting Fermi liquid at any fi- 
nite t should not be taken too far. In this regard, we 
should note that for \x\ 3> 2vt, the prefactor of the 
term oc (27ra;)~^ of the density correlation function (cf. 
Eq (21 )), which in equilibrium is proportional to the sys- 



tem compressibility, |57j remains equal to (minus) unity, 
which is the value that corresponds to a non-interacting 
Fermi gas (in an interacting Fermi liquid it would devi- 
ate from one). For t — -f-oo the same prefactor becomes 
(1 -f 7^) > 1, which does not relate easily to a non- Fermi 
liquid-like behavior. However, other correlation functions 
(like the one-particle density matrix discussed above) ex- 
hibit a similar behavior to the one-particle density ma- 
trix. For instance, let us consider the following correla- 
tors: 

C^{x,t) = (e2™</'(^.*)e-2™</'(0^*))^ (27) 

(28) 



CJ}{x,t) = (e™«(^^*)e-*"^(0'*) 



where (t){x) — ^ [<j)r{x) + <j)i{x)] and 9{x) = 
1 [(j>r{x) — (j)i{x)], the spatial derivatives of (p and 9 
are related to the (total) density and current density 
fiuctuations, respectively. [68^ Using exactly the same 
methods as above, we find (for L — >■ 00): 



c;'{^,tl 

{x) 



(0,m) 



c: 



C;{x,t) 
Cf-\x) 



RoV x^ 



2vt 



{2vtf 



Rq 

2vi 



2 2 
x'^ 



{2vtf 



m2(X2-l)/2 



(29) 
(30) 



The result of Eq. (26) is similar to the zero tem- 
perature result in equilibrium. However, the exponent 
of C^^{x,t — 00) is equal to 1 -I- 7^, and, even for 
an infinitesimal interaction {i.e. 7 ^ 1), it is always 
larger than the one that governs the ground state {i.e. 
equilibrium) correlations: [53l [55] 7^ = sinh^ 2/3(0) > 
7gq = 2sinh^/3(0). The reason for the larger expo- 
nent can be qualitatively understood from the following 



where C^^^^^^x) = AfJRo/x]^"''' and C^^'^'^x) = 
Af^lRo/x]"' /2 are the non- interacting correlation func- 
tion (where Af^ and A^^ are non- universal prcfactors). 
We note that the usual duality relation where (j> —?' 9 and 
K K~^, which one encounters when studying equilib- 
rium correlation functions, |57j still holds for these non- 
equilibrium correlators. Let us next analyze their asymp- 
totic properties. We consider only C2^{x,t), as identical 
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conclusions also apply to Cg{x, t) by virtue of the duality 
relation. For |a;| ^ 2vt, we have: 



2vt 



(31) 



Thus, up to the time-dependent pre-factor, correlations 
take the form of a non-interacting system of Fermions, 
C^^'"^\x). However, in the opposite limit (|x| ^ 2wi), 
this correlator exhibits a non-trivial power-law: 



Ro 



(32) 



Notice that this expression also describes the infinite- 
time behavior, which is controlled by an exponent equal 
to 771^(^2 + 1)^ being again different from the exponent 
exhibited by the same correlator in equilibrium, which 
equals 2m2 [cosh 2/3(0) - sinh 2/3(0)] = 2m^K. 

In order to understand why the behavior found for t — >■ 
oo in Eqs. ( 26|32 1 also holds for |x| <C 2vt, let us consider 
the initial state at zero temperature, po = |0)(0|. ITS] 
As mentioned above, this is a rather complicated ex- 
cited state of the Hamiltonian that performs the time- 
evolution, iff — Him- This means that, initially, there 
are a large number of excitations of Hi, namely, phonons 
with dispersion uj{q) — v{q)\q\. The distribution of the 
phonons {b^{q)b{q))o = sinh^/3(g) is time-independent 
and peaked at g = 0. Thus, within the approxima- 
tion where v{q) ~ v(0) = u, the excitations propagate 
between two given points with velocity v. Thus, if we 
consider the correlations at two points A and B sepa- 
rated by a distance the nature of the correlation at 
a give time t depends on wether the excitations found 
initially at, say, point A, have been able to reach point 
B or not. This is not the case if > 2vt, and thus 
correlations retain essentially the properties they had in 
the initial state. Thus, up to a time-dependent pref- 
actor, C^^{x,t) cx C^^\x). However, if the two points 
have been able to 'talk to each other' through the exci- 
tations present in the initial state, then correlations will 
be qualitatively different. This happens for a time t — to 
when the excitations propagating from A meet the exci- 
tations propagating from B, that is, for x — vIq = vto, or 
to = x/2v (we assume x > without loss of generality). 
Thus, for given separation x and time t, there is a length 
scale 2vt, which marks the transition between two differ- 
ent regimes in the correlations. In the instantaneous mo- 
mentum distribution, this reflects itself in a crossover as a 
function of time from a momentum distribution n{p) ex- 
hibiting a discontinuous Fermi liquid-like behavior, which 
is valid i.e. for \p\ <C {2vt)~^) to a power-law behavior of 

the form ~ Ipf^oP sgn(p), which applies for \p\ ^ {2vt)~^ 
but \p\ <^ R^^ (for \p\ ^ Rq^ we recover the free particle 
behavior corresponding to the Fermi-Dirac distribution 
function at T = 0). In the t — >■ cx) limit, by using the 
rcgularization scheme described above, the asymptotic 
momentum distribution at zero temperature can be ob- 
tained with the help of tables. [5^ The resulting formula 



behaves as the non-interacting Fermi-Dirac distribution 
for \p\ >• -Rq^, whereas for \p\ Rq^^ describes a 
non-Fermi liquid-like steady state: 



/(p,i ^ -foo) 



pRp 

2 



K,2_,i\pRo\)C-,2_,{\pRo\) 



+ K^{\pRo\)C^{\pRo\) , (33) 



where K^{z), C^{z) are the modified Bessel and Struve 
functions. [55] respectively. This expression yields a 
power law for |pi?o| ^ Ij where n(p,t — >■ -l-oo) ~ 
i — const. X {\pRo\)'' sign(p). Note that the momen- 



tum distribution n(p — 0, t) 



which is given by 



the invariance of the LM under particle-hole symmetry 

Let us finally present the generalization of the above re- 
sults for the one-particle density matrix to finite temper- 
atures, T > 0. For r < hvpRo (but T > Aq = 2nhv/L, 
so that we can neglect finite-size effects and effectively 
take the thermodynamic limit) C^^{x,t) takes the fol- 
lowing asymptotic form: 



C^^{x,t>0\T) = d°^{x\T) 



ttRq/X ^ 
dh{x\T) 

dh{x - 2vt\T)dh{x + 2vt\T) '^'^^ 
dh{2vt\T)dh{-2vt\T) 



(34) 



where C^°j(x|T) and dh{x\T) can be obtained from 

clp^{x\L) and d{x\L) by replacing L sm{Trx / L) /tt by 
A sinh(7ra;/A), where A = hvp/T is the thermal corre- 
lation length. At long times, hr{x,t\T) reduces to 



hr{x) 



ttRo/X 
sinh (ttx/A) 



(35) 



Therefore we again find that C^^{x, t — >■ oo|r) has a form 
similar to the the equilibrium correlation function at fi- 
nite temperature with a different exponent controlling 
the asymptotic exponential decay of correlations. No- 
tice that the exponential decay the correlations for i > 
is a direct consequence of the fact that the initial state 
has a characteristic correlation length, the thermal cor- 
relation length A — hvp/T.[69^ The exponential decay of 
correlations at finite T implies that the the steady state 
will be reached exponentially rapidly in a time of the or- 
der of h/T. It is also worth noting, however, that the 
above expression depends parametrically on the thermal 
correlation length A, and the Fermi velocity, vp, which 
enters in the expression for A = hvp/T, instead of the 
(renormalized) phonon velocity which enters in the ther- 
mal length Acq — hv/T, characterizing equilibrium cor- 
relations. Thus, since the velocity appears only through 
the definition of the thermal correlation length A, or, in 
other words, in combination with the temperature, the 
change from v to vp can be also understood as an change 
in the temperature scale. Furthermore, in a system with 
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Xp 

FIG. 2: Infinite time limit of tlie momentum distribution vs. 
Xp (A = hvp/T is tfie tfiermal correlation length in the initial 
state) for a non-interacting Luttinger model at finite temper- 
ature T that is quenched into an interacting state with re- 
pulsive interactions (corresponding to a Luttinger parameter 
K = 0.6). 




FIG. 3: Infinite time limit of the momentum distribution vs. 
Xp (X = hvp/T is the thermal correlation length in the initial 
state) for a non-interacting Luttinger model at finite temper- 
ature T that is quenched into an interacting state with at- 
tractive interactions (corresponding to a Luttinger parameter 
K = 1.7). 



Galilean symmetry, [5 7) we have that vK = vp and thus 
the parameter that controls the temperature scale now is 
the Luttinger parameter K, so that the asymptotic cor- 
relations at t oo can be regarded as the equilibrium 
correlations with a different exponent and an effective 
temperature, Tes = T/K. Thus, for repulsive interac- 
tions {i.e. K < 1) we could say that, besides modifying 
the exponent, the interaction quench increases the ef- 
fective temperature, whereas for attractive interactions 
(i.e. K > 1) the effective temperature is reduced after 
the quench. This effect has an impact on the momentum 
distribution at finite temperatures because it compen- 
sates the effect of the larger non-equilibrium exponent 
on the momentum distribution at finite T. To show this, 
we need to obtain the Fourier transform of hr{x). This 
can be done by relating it to an integral representation 
of the associated Legendre function P'^{z), [5S| and thus 
the Fourier transform of hr{x) can be written as: 



Kip) = 



X / nRr 



(7' + l)/2 



iXp 
'2ti 



X P] 



9 ~ 9 



2nRo 



A 



(36) 



Hence, the momentum distribution can be obtained by 
numerically evaluating the convolution with the Fermi- 
Dirac distribution function (cf. Eq. [Ts]) of the above ex- 
pression, Eq. (36 1. In Figs. [2] and |3| the momentum 



distribution of the interacting system for t ^ h/T is dis- 
played for a non-interacting LM that undergoes an inter- 
action quench with repulsive (corresponding to i^T = 0.6) 
and attractive (corresponding to K = 1.7) interactions, 
respectively. 



C. Suddenly turning-ofF the interactions 

Next we briefly consider the opposite situation to the 
one analyzed above, namely the case where interaction 
between the fermions in the initial state suddenly disap- 
pears. The fact that the initial state is a highly compli- 
cated state of the Hamiltonian that performs the time 
evolution (in this case H{ = Hq, cf. Eq. [T]), implies 
that we cannot expect that a Fermi liquid will emerge 
asymptotically at long times after the quench. Indeed, 
at zero temperature a thermodynamically large system 
approaches a steady state exhibiting equal-time corre- 
lations that decay algebraically in space. However, the 
exponents differ again from the (non-interacting) equi- 
librium ones. This can be illustrated by, e.g. computing 
the following correlation functions: 



\2vFt J 



m^{K-^-K)/2 



(37) 



= I^{x) (A) 



\2vFt J 

- (2vFif 



where 



X 

i?o 

X 



(38) 

(39) 
(40) 
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are the correlation functions in the initial (interacting) 
ground state {Ag/^ are non-universal prefactors). We 
note again that the duality 0^4) and K K^^ also 
holds in this case. The correlations in the stationary state 
that is asymptotically approached at long times read: 



lim CT{x,t) = Ae 



lim C^{x,t) = A4 



Ro 



X 

Ro 



X 



(41) 
(42) 



However, at short times, t ^ \x\/2vp, correlations look 
like those of in the initial state, up to a time-dependent 
prefactor: 



cr(.,^«N/2..) = (^) 



I^{x) (43) 



( 



(44) 



In this case the time-dependent prefactor has also a 
power-law form. 



III. LONG-TIME DYNAMICS AND THE 
GENERALIZED GIBBS ENSEMBLE 

Recently, Rigol and coworkers [TD] observed that, at 
least for observables like the momentum distribution 
or the ground state density, their long-time behavior 
following a quantum quench in an integrable system 
can described by adopting the maximum entropy (also 
called 'subjective') approach to Statistical Mechanics pi- 
oneered by Jaynes. [33 HEl IM] Within this approach, 
the equilibrium state of a system is described by a den- 
sity matrix that extremizes the von-Neumann entropy, 
S = — Trplnp, subject to the constraints provided by 
a certain set of integrals of motion of the system. In 
the case of an integrable system, if {/„i} is a set of cer- 
tain (but not all of the possible) independent integrals of 
motion of the system, this procedure leads to a 'general- 
ized' Gibbs ensemble, described by the following density 
matrix: 



1 

( 



(45) 



where Z^q ~ Tr ^^^^-^ . The values of the Lagrange 

multipliers A™ must be determined from the condition 
that 



{Im}gG — Tr [polm] — {In 



(46) 



where po describes the initial state of the system, and 
(•••)gG stands for the average taken over the general- 
ized Gibbs ensemble, Eq. (45 1. Although pi = |$(i = 



0))($(i = 0)1 in the case of a pure state, as was first used 
in Ref. \10i nothing prevent us from taking po to be an 
arbitrary mixed state and in particular a thermal state 
characterized by an absolute temperature T. In such a 
case, the Lagrange multipliers will depend on T or any 
other parameter that defines the initial state. 

Rigol and coworkers tested numerically the above con- 
jecture by studying the quench dynamics of a ID lattice 
gas of hard-core bosons (see Ref.[Tni[ISlfor more details). 
The question that naturally arises then is whether corre- 
lations of the LM also relax in agreement with the above 
conjecture. In other words, does the average {0){t) at 
long times relax to the value (O)gG = Tr pgQ O, for any 
of the correlation functions considered previously? We 
shall first discuss an observable for which the generalized 
Gibbs ensemble fails to reproduce their expectation val- 
ues. Moreover, by considering the correlation function of 
the current operators {i.e. O = Jr{x), cf. Eq 22 1, we will 



illustrate why it works in the LM. Calculations of other 
observables can be found in the appendices. 

Let us define the generalized Gibbs ensemble for the 
Luttinger model (LM). Since the final Hamiltonian (in 
the TV = and J = sector) is diagonal in the b bo- 
son operator basis, i.e. Hlm = X^g^o ^Hq)Hi)^ 
a natural choice for the set of integrals of motion is 
Ijn — I{q) = n{q) = {q)h{q) for all g 7^ (a more 
complete version of the ensemble should also include N 
and J, but this will not be necessary as we work in the 
thermodynamic limit). Thus, for the quench from the 
non- interacting to the interacting state (cf. Sect. II B I, 
where the initial state is |$(i — 0)) = |0), the Lagrange 



multipliers are determined by Eq. (46), which yields: 



{i{q)). 



{'n{q))gG = sinh^ 13 (q) 



1 



(47) 



Indeed, this result can be quickly established by realiz- 
ing that pgG has the same form as the density matrix of 
a peculiar canonical ensemble where the temperature on 
each eigenmode of the final Hamiltonian depends on the 
wave-vector q, that is, T{q) = hv{q)\q\/ X{q). Alterna- 
tively, one can also regard it as an ensemble where the 
effective Hamiltonian that defines the Boltzmann weight 
is given by Hes/Tes = J2q^Q Hq)n'{q)- However, it is 
worth noting that is diagonal in n{q), and therefore 
it does not capture the correlations existing in the initial 
state between the q and —q modes. Mathematically, 

{n{q)n{-q)) = sinh^ /3(g) cosh 2/3(g) {n{q)n{-q))gG 
= {n{q))gG{n{-q))sG - sinh^ ^{q) . (48) 

As matter of fact, since n{q)n{—q) commutes with 
we conclude from the above that {n{q)n{—q)) does not 
relax to the value predicted by Pgc- Although this de- 
fect of pgG can be fixed by enlarging the set of integrals 
of motion to include the set I'{q) = n{q)n{—q) as well, 
and the corresponding lagrange multipliers, we shall show 
below that this is not needed. The reason is that the 
correlations missing in the generalized Gibbs ensemble 
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as defined above yield a vanishing contribution in the 
thermodynamic limit to the simplest correlation func- 
tions. However, before discussing this point, it is worth 
mentioning one important exception to the class of ob- 
servables that relax according to the generalized Gibbs 
ensemble, namely the squared fluctuations of the energy: 



a2 = (ij2) „ 



P,<1 



huj{p)huj{q) 



which yields 



X [{n{p)n{q))-{nip)){n{q))] (49) 
= 2<^Ig = E,sinh2 2/3(g);i2^(<j)2. 



Again, since the operator is conserved, violates the 
relaxation hypothesis. However, it is tempting to argue 
since (like (H)) is a non-universal property of the LM 
model, this violation is less problematic than a violation 
in the asymptotic behavior of the correlation functions 
would be, as the latter tends to be more universal. 

In order to proceed with the discussion of the validity 
of the generalized Gibbs ensemble, let us consider, for 
the sake of simplicity, the correlations the current opera- 
tor. In particular, we shall study the following two-time 
correlation function (no time ordering is implied): 

Cj,ix,t,T) = {Jrix,t + T/2)JriO,t-T/2))^. (50) 

where (. . .)t stands for average over the initial thermal 
ensemble described by pi — e~^'/"^/Zo, with Hi = Hq 
(cf. Eq.[T]). Using ([l5]), and ([l6]), we obtain: 

(27r)2 L 

^ ' q>0 ^ ' 



X 



hT/2)r(9,t-r/2) [l + nB{q,T)] 
T/2)g{q,t-T/2)nB{q,T) 



+e-^'^y*{q,t + r/2)f{q,t- 
+e-^'^^giq,t + T/2)g*{q,t~ 



T/2)nB(g,r) 

t/2) [l + nB{q,T)]}, 



(51) 



the initial Bose distribution of modes 



being nB{q,T) 

(^0(9) = VF\q\) 
and where / and g are defined in Eqs. ( ASP and ( A9 1. In 
the following we shall argue that, in the limit t — > +00 
the above expression reduces to the following correlator 
in the generalized Gibbs ensemble: 



cf (x, t) ^ Tr [pgG(r) Jr{x, t)MO, 0)] 



(52) 



where pgQ,{T) is the extension to an initial thermal state 
of the generalized Gibbs ensemble introduced above (no- 
tice that since [Hi, I{q)] — and therefore [Hi, pgc] — 0, 
it is in principle possible to define time-dependent corre- 
lation functions on this ensemble, just as they are defined 
in the equilibrium ensembles). For this (thermal) initial 
condition (461 fixes the values of X{q,T), which now de- 
pend on f3{q) and on the temperature T: 



sinh" /3((7) [1 + UBiq, T)] + cosh' /3(g)nB(g, T) 

_ 1 



(53) 



Introducing this result into the mode expansion for 
Eq. ( 52 ) , we arrive at 



^ ^ q>0 ^ 



-qao 



X {e*9(^-"^) cosh' f3{q) [1 + {n{q))] 
+ e*«(^+^^) sinh' I3iq) {n{q)) 
+ e-*«(^~^^) cosh' ^(q) (n(g))} 
+ e-*«(^+^^) sinh' /3(q) [1 + {n{q))] . 



(54) 



Using the expressions for f{q,t) and g{q,t) given in the 
Appendix, it can be shown that, in the limit t — ^ -\-oo 
of Eq. (51), the rapidly oscillating terms that depend 
only on t can be dropped, and therefore Eqs. (51) and 



(54) become the same. Thus, the above current correla- 



tion function relaxes according to the generalized Gibbs 
ensemble. 

Let us close this section with a digression of higher 
order current correlation functions. In principle, these 
correlations depend on correlations between n(q) and 
n{—q) {i.e. on {n[q)n{—q))), which exist in the initial 
state. These are not captured by the above simple- 
minded generalized Gibbs distribution, which only con- 
tains information about the expectation value of n(g) in 
the initial state. However, a simple argument based on 
momentum conservation and counting powers of L (the 
system size) shows that the contribution of these cor- 
relations vanishes in the thermodynamic limit, as men- 
tioned above. To illustrate this point, let us consider 
computing the four point current correlation function, 

{JriXi,ti)Jr{x2,t2)Jr{x3,t3)Jr{x2,t4))T. UpOn USing 

the mode expansion for the current operator Jr{x,t) — 
d(f>rix,t)/2Tr, the correlation function can be expressed 
as a sum over expectation values of the mode creation 
and destruction operators, the b operators. However, it 
is important that, since each of these mode expansions of 
the current operator carries a factor of , the four point 

correlation function is thus proportional to L^' . Upon 
applying the standard Wick's theorem to the expectation 
value of the b operators, momentum conservation requires 
that at least two of the four momenta being summed over 
must be equal. In general, the two independent momenta 
are not equal, and therefore we get a finite contribution in 
the limit L — > cx), which involves only a product of expec- 
tation values of the operator n{q). However, when two 
independent momenta being summed over coincide, the 
expectation value becomes {n{q)n{—q))T and thus, in the 
thermodynamic limit, the sum turns out to be of order 
j^, therefore vanishing as L — >■ 00. This argument can be 
extended to higher order correlations of the current op- 
erator. In a sense, it also justifies the use of Wick's the- 
orem when computing higher order correlations (in the 
thermodynamic limit) using the above 'simple-minded' 
generalized Gibbs ensemble. 
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IV. RELEVANCE TO EXPERIMENTS 

As we described in the introduction, ultracold atom 
systems are the ideal arena to study quench dynamics. 
This is because they are, to a good approximation, iso- 
lated systems. Furthermore, as far as one dimensional 
systems are concerned, there are already a number of 
experimental realizations, including experiments where 
quench dynamics has been already studied. [13J [501 EI] 
Thus, in this section we would like to discuss the possible 
experimental relevance of the results obtained in previ- 
ous sections. As mentioned above, this must be done with 
great care because our results have been obtained using 
a field theory model, namely the Luttinger model (LM), 
which can be regarded as 'caricature' of the Hamiltoni- 
ans describing real systems of ultracold atoms confined 
to one dimension. We must emphasize that the situa- 
tion in the case of quantum quenches in particular, and 
of non-equilibrium dynamics in general, is very different 
from the analysis of low-temperature phenomena in equi- 
librium. In the latter case, the experimental relevance 
of models such las the LM is well established by using 
renormalization-group arguments. This has been put to 
test over the years using a large variety of numerical and 
also (when possible) analytical methods. By contrast, in 
the case of non-equilibrium dynamics we travel through a 
largely uncharted land, and much needs to be understood 
in order to achieve a similar level of rigor as in the equilib- 
rium case. Thus, it is convenient to regard models as the 
LM as 'toys', which can provide us valuable lessons and 
insights into non-equilibrium dynamics of strongly cor- 
related systems. With this cautionary remarks, we can 
proceed to discuss some experimental systems for which 
the above results could be of some relevance. 

As mentioned in Sect. [lljthe LM is the exactly solv- 
able model describing the renormalization-group fixed 
point of a general class of interacting one-dimensional 
models'^, known as Tomonaga-Luttinger liquids. This 
class includes systems such as the one-dimensional Bose 
gas interacting via a Dirac-delta potential (which is solv- 
able via the Bethe-ansatz [J?]) as well as many other 
systems of interacting Bose gases with repulsive inter- 
actions (such as dipolar) or Fermi gases with both at- 
tractive and repulsive interactions. With the caveats of 
the previous section, it would be interesting to test the 
results obtained using the LM in one of these systems. 
However, the dynamics may be strongly modified by the 
fact that higher energy states will be also excited fol- 
lowing a quantum quench. Such higher energy states are 
not, in general, correctly described by the LM. The worse 
situation may be expected in the case of a system where 
interactions are short ranged, because, at least from a 
perturbative point of view, an interaction quench will 
couple excited states with the same strength. Thus, one 
possible way around this problem when trying to com- 
pare with results for the LM model would be to study 
experimentally (or numerically) interaction quenches in 
systems with long-range interactions. This system pro- 



vides us with a much more faithful realization of the LM, 
which involves long-range interactions. Furthermore, the 
sudden connection of interactions is likely not to scat- 
ter particles to high energy states because interaction 
decreases with the momentum exchanged between the 
colliding particles. One system satisfying this require- 
ment is a single-species dipolar ID Fermi gas confined to 
one dimension by a strongly anisotropic trapping poten- 
tial. Since for a single species Fermi gas contact inter- 
actions in the p-wave channel are weak away from a p- 
wave Fcshbach resonance, the dominant interaction is the 
long-ranged dipolar interaction when the dipoles are all 
aligned by an external field (electric or magnetic, depend- 
ing on whether the dipole is electric, like in hetero-nuclear 
molecules, or magnetic, like in Chromium). When con- 
fined to a one-dimensional trap, the dipolar interaction 
between the atoms can be approximated by the potential: 



1 



47reo (a;2 i?2)3/2 



(55) 



where D is the dipolar momentum of the atoms, 9 is the 
angle subtended by the direction of the atomic motion 
and the polarizing field, and X{9) = (1 — 3cos0). Since 
in this case g2('?) = 34(9) X{9), a sudden change in 
the interactions can be produced by a sudden change in 
alignment of the field with the direction of motion, that 
is, a change in 9. In particular, a change in 9 away from 
the value 9m — cos~^(^) would lead to suddenly switch- 
ing on the interactions between the fermions. [3] At zero 
temperature, the momentum distribution f{p,t) (which 
can be probed by time of flight measurements) following 
the quench into the interacting system would evolve as 
described in Sect. II B (cf. Fig. [I]), with the discontinuity 
at the Fermi level dying out as t~''^ . However, currently 
atomic gases are produced at temperatures T ~ 10% 
to 20% of the Fermi energy, and this would complicate 
the observation of this effect. If much lower tempera- 
tures could be reached in experiments, so that the ap- 
plication of the LM becomes much more reliable, we ex- 
pect that in a time of the order of h/T the quenched 
dipolar gas reaches a stationary state characterized by 
a momentum distribution that differs from the thermal 
one. However, the calculations of f(p, t) presented in 
Sect. II B| (cf. Figs. [2] and |3| show that the differences 
between the non-equilibrium and equilibrium results in 
the stationary state may be well below the current ex- 
perimental resolution. Alternatively, instead of measur- 
ing the momentum distribution, one can try to deter- 
mine the non-equilibrium exponents by measuring noise 
correlations in the time-of-flight images |62j or through 
interferometry. |63j 

To sum up, we have studied correlations in the Lut- 
tinger model (LM) following an interaction quench. By 
studying the situations where the interactions are sud- 
denly turned on and off, we have shown that the correla- 
tion dynamics is dominated by the initial state correla- 
tions at short times. However, in the thermodynamic 
limit, the LM reaches a non-equilibrium steady state 
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characterized by a set of non-equilibrium exponents that 
differ from the exponents that describe the decay of the 
same correlations in equilibrium. This behavior can be 
obtained from a 'simple-minded' generalized Gibbs en- 
semble, [ini US US] The correlations also exhibit a 
'light-cone effect', which was previously observed in cor- 
relations following a quantum quench from a gapped to 
a critical state. |T2l[TS]. Finally, we have discussed the 
relevance of our results for the LM for experiments with 
ultracold atomic gases confined in highly elongated (one- 
dimensional) traps. We have argued that the most appro- 
priate scenario for the observation of the effects described 
here may be a one-dimensional dipolar Fermi gas. 



evolution of b{q) for t > 0. For Hamiltonians like (Al) 



this can be done exactly because H{ = H{t > 0) can 
be diagonalized by means of the canonical Bogoliubov 
( "squeezing" ) transformation: 

a{q) = cosh/3(g)6(g) +sinh/3(g)6l'(-9). (A3) 

Upon choosing 

.9(9) 



tanh 213 (q) 



(A4) 



ujoiq) + nT'{q)' 
the Hamiltonian at t > is rendered diagonal: 

H, = H = Eo + J2f^^il)^HqHq), (A5) 
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Appendix A: Quadratic Hamiltonians 

In this appendix we study the case of a quantum 
quench in a model described by a quadratic Hamiltonian: 

H{t) = Y,h [c.o(9) + m(<z, t)] 6t {q)h{q) 

9 

+ \Y.^g{<l,t) [b{q)b{-q)+hHq)b\-q)] , (Al) 
9 

where [&(q), 6^(q')] = (5g,g', commuting otherwise. We 
will assume that the quench takes place at t = 0, so that, 
within the sudden approximation, the system is described 
by Hi = H{t < 0) for t < and by Hf = H{t > 0) for 
t > 0. Furthermore, in order to simplify the analysis, 
we assume that m{q,t ^ 0) = g{q,t ^ 0) = 0, and 
m{q,t > 0) = m{q) and g{q,t) — g{q). Notice that the 
initial Hamiltonian is diagonal in the b operators: 



Hi^Ho = Y^hu:o{q)b\q)b{q). 



(A2) 



In order to obtain the time evolution of operators O = 
0[{b'^[q),b{q)}] after the quench, we recall that, in the 
Heisenberg picture, 0[t > 0) = ^iHit/hQ^-iHft/h ^ 
0{{b{q,t),b'^ {q,t)}), and therefore all that is needed to 
solve the above quench problem is to obtain the time 



where Eq is the energy of the ground state of H (relative 
to the ground state energy of Hg) and 



ojiq) = ^[ojoiq) +m{qf~[g{q)]^ (A6) 

the dispersion of the excitations about the ground state 
of H{. The evolution of the a{q) is given by a{q, t) = 
^iHtt/h^^^q-^^-tHit/h ^ e-''^(9)*a(g). By application of a 
direct and reverse Bogoliubov transformation, one can 
obtain the time evolution of &((?): 

b{q, t) = f{q, t) b{q) + g* {q, t) b^-q), (A7) 



where 



f{q,t) = cosuj{q)t — ismoj{q)t cosh2/3(g), (A8) 
g{q,t) = ismuj{q)t sinh2l3{q). (A9) 



It is easy to check that (A7) obeys the initial condition. 



b{q,t = 0) = b{q), and also respects the equal-time com- 
mutation rules, 

[b{q,t),b{q',t)]^{f{q,t)g*{q,t) 

-g*{q,t)f{q,t))S,^^,^0, (AlO) 

Hq,t),b^{q\t)]^{\f{q,t)\'-\giq,t)\') S,,,, 

= S,.^,' (All) 

Thus, a quantum quench described by a quadratic Hamil- 
tonian can be solved by means of a time-dependent 
canonical transformation. 

When the quench is reversed, i.e. when the case with 
m{q,t ^ 0) = g{q,t > 0) = 0, and m{q,t < 0) = m{q) 
and g{q, t < 0) = g{q) is considered, the roles played by 
the initial and final Hamiltonians are also reversed: the 
final Hamiltonian is now diagonal in the &'s, H{ = Hq, 
whereas the transformation of Eq. ( A3 1 renders diagonal 
the initial Hamiltonian, Hi = H. Therefore, in this case 
the evolution of the b operators is trivial: Hf -. b{q, t) — 
g-«i^o(9)*^ whereas the evolution of the a's is given by 



a{q,t) ^ fo{q,t)a{q) + gl{q,t) aH~q), 



(A12) 



where 



fo{q,t) ~ cos ujo{q)t — ismujo{q)t cosh2(3{q), (A13) 
go{q,t) = -ismu!o{q)t sinh2/3(q). (A14) 
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Appendix B: Details of the calculation of the 
one-particle density matrix in the Luttinger model 

In this Appendix, we shall provide the details of the 
calculation of non-equilibrium one-particle density ma- 
trix: 



(Bl) 



To this end, the formula ( 14 1 is used. In normal ordered 
form: 



/L 



(B2) 



where the normal order is defined: 



The boson field ^a{x) is given by Eq. (16). Hence 



X : e 



-i[<Pr(x)~4>Ao)] 



, (B4) 



where we have used the identity e^e^ = el"^'^' e^e^, 
which holds provided [A, B] is a c-number. Using that 
(ao — )■ 0"*" is the short-distance cut-off): 



[<I>,(a:),<ft(0)]=5^(^|e -e^ 

q>0 ^'^ ' 

_ _ \^ _ g-2irao/-Lg2iir£c/L 

we arrive at the following expression for C^^, (a;,i): 
where 



1 



1L sin[f (x + ia)] 
Fr{x, t) = e'"'/'' [<S>r{x) - $r(0)] e"'"''/" 

g>0 ^ 



(B5) 
(B6) 

(B7) 

(B8) 
(B9) 
(BIO) 



X [f{q,t)biq)+g*{q,tp{-q)] . 

(Bll) 



To derive the last expression we have used Eq. ( |A7| . 
Employing the identities e^e^ = gl^'B]/2^A+B (p^QyiJed 

[A,B] is a c-number) and that {e^) = e<-°> + 5<(^-<-°»'> , 
we obtain: 



(B12) 



where we have used that {[F} {x , t)]'^) = {F^{x,t)) = 
because ([6^(9)]^) = {[b{q)]^) = {b^-q)b{q)) = 0, and 
since the commutator [F^{x, t), Fr{x, t)] is a c-number, it 



can be safely replaced by {[F^{x,t),Fr{x,t)]). Note that 
for t = Fr{x,t = 0) contains only b{q) and thus the 
average {Ff{x, t)Fr{x, i)) = at T = 0. In Eq. ( |B12| ) the 
exponent can be expanded to yield: 



g>0 



{Fl{x,t)Fr{x,t))^Y.['^]^ 



qL 



-qao\„iqx _ 1 |2 



X [\f{q,t)\^nB{q) + \9{q,t)\^{nB{q) + l)], (B13) 

being nsiq) = {a^{q)a{q)) = (e"^l«l - 1)"^ the distri- 
bution of Tomonaga bosons in the initial state (which 
has been assumed to be a mixed thermal state), and 
A = hvp/T is the thermal correlation length. We next 
evaluate explicitly the above result in several limiting 



1. Zero temperature and finite length 

Let us now consider the T = limit of the above cor- 



relator, where nsiq) = 0, and thus, using (A8I and (A9|, 
Eq. (B13) simplifies to: 



{F^{x,t)F,{x,t)}T=o = (^) sinh2[2/3(g)] 

q>0 ' 

X (1 - cosgx) [1 - cos2t;(g)|g|t]. (B14) 

To make further progress, we assume that sinh2/3(g) ~ 
^g-|9flol/2 -vifhere i?o is the range of the interaction. Fur- 
thermore, we replace viiq) by w = ?;(0) what allows 
us to safely take the limit oo — > 0+. Next, in order to 
simplify the computation, we introduce the quantity 



(?>0 



£^{z) = ^ ^ ^ ) e-«^« cosgz. 



(B15) 



which can be readily computed to give 

7ri?o 



In 



L 



d{z + iRo\L) 



L 



In 2, (B16) 



where d(z\L) — L \ sin(7rz/L)|/7r is the cord function. Us- 
ing this result into Eq. (B14), yields the following ex- 



pression for the one-particle density matrix: 

d{iRo\L) 



C,p^ix,t>0\L)^Gi'>\x\L) 



d{x + iRo\L) 
d{x + 2vt + iRQ\L)d{x - 2vt + iRo\L) '^'^^ 



d{2vt + iRo\L)d{-2vt + iRo\L) 



(BIT) 



Taking into account that Rq/L <^ 1 and, we obtain the 
result quoted in the main text, Eq. (20), in the scaling 
limit. 
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2. Thermodynamic limit and finite temperature 



We note that the result of Eq. ( B26 1 can be obtained from 



We next consider Eq.(B12) for i — > oo and finite tem- 
perature, T. Equation ( B13| ) can be recast as: 



{Fl{x,t)Fr{x,t))T = {Fr{x,t)Fr{x,t))T=0 

+ 'H{x)+g{x,t) (B18) 
where we have introduced the following functions: 

2 I [1- cos qx)nB{q), (B19) 

Jo 1 

°° dq 



nix) 

g{x,t) = 2j^ f 

Jo 



-qRo 



(1 — cosqx) 



[1 " cos{2vqt)] nsiq), 



(B20) 
(B21) 



which hold in the thermodynamic limit and upon replac- 
ing v{q) by i; = v{q = 0) and sinh2/3(g) = -ye'^'^^''^/'^ 
as we did in the previous section. We next define the 
function 



Jo <1 e^P - I 

which can be evaluated to yield: [58] 

r(i + A-V) 



g{u; r) = 2 In 



r[l + A-i(r + m)] 



(B22) 



(B23) 



where T{z) is the Gamma function. In the limit where 
r <C u, and using that|58j r(z)r(l — z) ~ 7r/sin(7rz), the 
above expression reduces to 



l{u; r) = - In 



dh{ir\T) 
dh{u + ir\T) 



-In 


u + ir 




r 



In the previous expression we have defined: 



A 



dh{z\T) = — |sinh (ttA """z) | 



(B24) 



(B25) 



Combining this result with Eqs. ( |B18| )-( |B2l| andjB12|, 
it is seen that the second term in Eq. (B24) ex- 
actly cancels the contributions from G^"-' (x\L) and 
{Fr{x,t)Fr{x,t))T=o in the thermodynamic limit, and 
therefore, 



C^Jx,t\T) = Gi'^HxlT) 



dh{iRa\T) 



dh{x + iRo\T)_ 
dh{x + 2vt + iRo\T)dh{x - 2vt + ii?o|T)"l '''^^ 



dh{2vt + iRo\T)dh{-2vt + iRo\T) 



(B26) 



where 



G(°)(x|T) = 



i ttA-i 
27r sinh [vrA^^ (x + iao)] 



Eq. (B17) upon making the replacement Lsin(7rL~^a;)/7r 
by Asinh(7rA^^a;). Taking into account that Rq/L ^ 1, 
we retrieve the result quoted in the main text, Eq. p4[ ), 
in the scaling limit. 

Appendix C: One-body density matrix of the 
Luttinger model in the generalized Gibbs ensemble 

Next we take up the calculation of the one-body den- 
sity matrix in the generalized Gibbs ensemble for the 
Luttinger model discussed in Sect. |III| That is, we shall 
evaluate the expression at T = 0. 



(CI) 



Using the bosonization identity, Eq. ( |14[ ), we can write 
the expression as follows: 



Taking into account that 



Frix) = <i>r{x) - $r(0) 



E 

g>0 



27r 
qL 



1/2 



(C2) 
(C3) 

(C4) 
(C5) 



X [cosh f3{q)a{q) - sinh/3(q)a^(-g)] . (C6) 

The expression for C^{x) can be easily computed by 
using the trick of regarding pgc as a canonical ensemble 
with g-dependent temperature. Thus, following the same 
steps as in the previous section we arrive at: 



Given that 

{F^{x)Frix))gG = sinh^ 2/3 [2?,(0) - Vr{x)] 
where 

'2^~ 



Vr{x) =Re<^ 
U>o 



qL 



In sin — {x + iRo) 



In 2 



L 



Hence, taking the thermodynamic limit 

i Ro 



(C7) 
(C8) 

(C9) 
(CIO) 

(Cll) 



27r(x -|- ia) 

Thus we see that one recovers the same results as 



(B27) Mnvt^+oo C^Ax,t), Eq. (23) 
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146) to Statistical Mechanics. See for instance the critique 
by Ma in Ref. 1641 However, the maximum entropy ap- 
proach is advocated by Balian. ^\ 
[67] An exception are pathological cases, like the Coulomb in- 
teraction, where both v{q)/\q\ and sinh2/3(g) are singular 
for q = 0. 

[68] In a Luttinger liquid, the correlator describes the 
fluctuations of wave number close to 2mfc_F (where kp 
is the Fermi momentum) of the density-correlation func- 
tion |55l 157) . In the presence of a periodic potential of 
periodicity equal to 2mkF the system may become an 
insulator. [57] The power-law behavior exhibited at zero 
temperature, and in the thermodynamic limit) is a con- 
sequence of the gapless spectrum and the absence of long 
range order in the density. In the insulating (i.e. gapped) 
phase, this correlation function decays to a non-zero con- 
stant at long distances, which is a consequence of the exis- 



tence of long-ranged order in the density at wave number 
2mkF- Similarly, in equilibrium Cg measures the phase 
fluctuations, and exhibits a power-law, reflecting the ab- 
sence of long range order in the phase. However, in the 
Luttinger model of interest for us here, terms with m > 
are absent from the density operator, which is given by 

p{x) = Jr(x) + Jl{x). [53 

[69] A similar situation is found when analyzing quenches at 
T — from a non-critical (that is, gapped) into a critical 
(that is, gapless) state. In that case, the role of A will 
be played by the correlation length of the system that 
is determined by the (inverse of the) energy gap in the 
initial state. |T2l |48] 

[70] We assume implicitly that both j3{q) and v{q) are not 
singular at g = 0. 



